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ABSTRACT 

New spectroscopic observations show that the double degenerate system NLTT 16249 is in a close 
orbit (a = 5.6±0.3 Rq) with a period of 1.17 d. The total mass of the system is estimated between 1.47 
and 2.04 M Q but it is not expected to merge within a Hubble time-scale (t mcrgc w 10 11 yr). Vennes 
& Kawka (2012, ApJ, 745, L12) originally identified the system because of the peculiar composite 
hydrogen (DA class) and molecular (C2-DQ class-and CN) spectra and the new observations establish 
this system as the first DA plus DQ close double degenerate. Also, the DQ component was the first of 
its class to show nitrogen dredged-up from the core in its atmosphere. The star may be viewed as the 
first known DQ descendant of the born-again PG1159 stars. Alternatively, the presence of nitrogen 
may be the result of past interactions and truncated evolution in a close binary system. 
Subject headings: binaries: close - stars: fundamental parameters - white dwarfs 



1. INTRODUCTION 

The high-proper motion star NLTT 16249 is a double 
degenerate system showing hydrogen lines and molec- 
ular c arbon and cyanogen ba nds in its optical spec- 
trum (|Vennes fc Kawkal 120121 ). Al though large ra- 
dial v elocity variations were noted by iVennes fc Kawkal 
(2012), the orbital parameters, period and separation, 
are yet to be determined. The detection of photo- 
spheric nitrogen in the carbon-rich (DQ) component 
of this system was a first occurrence for this class of 
objects and, according to the dred ge-up scenario com- 
monly applied to DQ white dwarfs jPelletier et al.lll986l : 
iMacDonald. Hernanz. fe JoseHl99"8l ). it implied the pres- 
ence of nitrogen in the white dwarf core. The other com- 
ponent is a hydrogen-rich (DA) white dwarf with a mass 
above average but with a luminosity similar to that of 
the DQ white dwarf. 

The two characteristics of the system, a nitrogen- 
enriched DQ component and a likely close orb it, may 
or may not be related. Alt haus et al.1 (|2005[ ) estab- 
lished a clear evolutionary link between born-again 
stars and a nitrogen/oxygen enrichment in PG1159 
stars and their DQ descendents. Members of the 
PG1159 class have a helium-rich surface with no- 
table enrichment in carbon, nitrogen and oxygen 
dWerner. Heber. fc Hunger]ll991l:lDreizler fe Heberl ll998l: 
I Werner fe Herwigl I2006D . I Althaus et al.l (|2005r )argued 
that not all DQ stars followed the born-again path 
because the carbon abundance is lower in some of 
these objects than would be expected following this 
evolutionary path. H owever, the logC/He versus 
T R ff tr end established bv | Dufour. Bergeron, fe Fontaine! 
(2005) and lKoester fc Knist[ (2006') does a ppear to follow 
the model predictions of lAlthaus "eTaLl (l2005lh More- 
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over, iDufour. Bergeron, fc Fontaine! (|2005l) noted two 
discernible tracks with one at a higher carbon abundance 
that was attributed to a thinning of the helium layer. 
Therefore, data and born-again models may be compat- 
ible for the bulk of DQ white dwarfs. However, nitrogen 
is not detected in any DQ stars other than in the DQ in 
NLTT 16249. 

In this Letter, we present convincing evidence that 
the components of the binary NLTT 16249 are in a 
close orbit. Our new radial velocity measurements (Sec- 
tion 2) and our origin al model atmosphere analysis 
(|Vennes fc Kaw ka 2012) help constrain the component 
properties and offer clues to the origin and evolutionary 
prospect of the system (Section 3). We summarize and 
discuss these new results in Section 4. 

2. OBSERVATIONS 

IVennes fc Kawkal (|2012f ) described the discovery 
and initial observations of the double degenerate 
NLTT 16249: on UT 2010 March 26, they obtained three 
low-dispersion spectra using the R-C spectrograph at- 
tached to the 4 m telescope at Kitt Peak National Ob- 
servatory (KPNO). Next, on UT 2010 November 7, they 
obtained an echelle spectrum using the X-shooter spec- 
trograph attached to the UT2 at Paranal Observatory. 

The apparent radial velocity variations and the com- 
posite nature of the original spectra prompted us to 
obtain a more comprehensive time series. We re- 
observed NLTT 16249 at Siding Spring Observatory 
(SSO), KPNO, and the MDM Observatory (on Kitt 
Peak). 

On UT 2011 December 2-3, we obtained a set of spectra 
with the Wide Field Spectrograph (WiFeS, iDopita et al.l 
I2007T ) attached to the 2.3m telescope at SSO. On De- 
cember 2, we obtained three exposures of 15 minutes 
each, and, on December 3, two exposures of 30 min- 
utes each. We used the B3000 and R7000 gratings which 
provided spectral ranges of 3200-5900 A at a resolution 
R = A/AA = 3000, and 5300-7000 A at R = 7000, re- 
spectively. We used the RT560 dichroic beam-splitter to 
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Radial velocities. 
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Fig. 1. — Frequency analysis (inverse of the reduced Xr versus 
frequency, l/P) of the Ha data (i>da> Table[T) including the KPNO 
data only (long-dashed line), the KPNO and SSO data (short- 
dashed line), and the complete data set (X-shooter, KPNO, and 
SSO; full line). The horizontal lines mark lcr confidence levels for 
the corresponding data sets. 

separate the incoming light into its red and blue com- 
ponents. We maximized the signal-to-noise ratio of each 
observation by extracting the spectrum from the most 
significant (< 6) traces. Each trace was wavelength and 
flux calibrated prior to co-addition. The spectra were 
wavelength calibrated using NeAr arc spectra that were 
obtained following each observation. 

Next, on 2012 January 4-6, we used the 4 m telescope 
and R-C spectrograph at KPNO. The KPC-24 grating 
in second order and the T2KA CCD delivered a spec- 
tral resolution of ~0.9A from 6030 to 6720A with the 
slit width set to 1.5". We sorted the first and second 
orders using a GG495 filter. We obtained 17 exposures 
of 30 minutes, each one followed by a comparison arc 
(HeNeAr). 

Finally, on 2012 January 18, we used the MDM 2.4 m 
Hiltner telescope and modular spectrograph (modspec). 
A thick 2048 2 CCD gave 1.7 A pixel" 1 and ~ 3.5 A res- 
olution from 4460 to 7770 A, with vignetting toward the 
ends of the range. We obtained three consecutive expo- 
sures (5, 5, and 10 minutes). All spectra were reduced 
using standard IRAF procedures. 



3. ANALYSIS 
3.1. Radial velocities and period analysis 

The narrow Ha line core proved suitable for precise ra- 
dial velocity measurements of the DA white dwarf. When 
available, the sharp edges of the C2 and CN band heads 
also offer means to measure relative DQ velocities. Con- 
sequently, we measured the DA radial velocity by fitting 
Gaussian functions to the narrow Ha core (uda in Ta- 
ble [T]) . The measurement errors were estimated by vary- 
ing the location of the continuum when fitting the line 
core and we noted fluctuations of 2-5 kms" 1 depending 
on the signal-to-noise ratio. The centering of the tar- 
get on the slit is the largest source of systematic errors. 
These errors were estimated at ~ 1/10 of a resolution el- 
ement or ~ 5 kms" 1 for a resolution R k, 7000. On the 
other hand, residuals in the wavelength calibration never 
exceeded 1-2 kms -1 and were similar to velocity fluctua- 
tions of the O IAA6300.23, 6363.88 sky lines. Adding sys- 
tematic and random errors, we conservatively estimated 
the total error at 10 kms -1 . 

Next, we re-measured the DQ velocity in the UVB 
X-shooter spectrum using the C2 Swan band head at 
5162. 2A and the CN violet band head at 3883.4A and 
listed the average measurement in Table Q] The mea- 
su rements are consistent w ith those originally employed 
bv lVennes fc Kawkal (|2012f ). The MDM and SSO spectra 
were then cross-correlated with the X-shooter spectrum 
to extract new DQ velocities (wdq in Tabled]). Consecu- 
tive spectroscopic series obtained at SSO or MDM were 
co-added to maximise the signal-to-noise ratio and de- 
crease the error in veloc ity measurements . The original 
KPNO data obtained bv lVennes & Kawkal (|2012h are not 
included in the present set because of their lower spectral 
resolution (R « 1000). 

Finally, we fitted the sinusoidal function v — 7 + 
K x sin (27r(£ - T )/P) as a function of time (HJD t) 
to the DA radial velocity data and determined simul- 
taneously the initial epoch (To), period (P), mean ve- 
locity (7) and velocity semi-amplitude (K). We nor- 
malized the x 2 function so that the minimum reduced 
Xrmin = The la range is then estimated using 
X 2 r,ia = Xlmin + c/{N - p) = 1 + 4.7 /(N - 4), where 
N is the number of measurements, p = 4 is the number 
of parameters, and c = 4.7 is the appropriate constant for 
a la error range and 4 parameters. Figure [T] shows incre- 
mental frequency (l/P) analyses. The period extracted 
from the KPNO radial velocity data alone (N = 17), 

P= 1.161 ±0.026d, 

is sufficiently accurate to constrain the DA velocity am- 
plitude and mass function of the DQ companion. Adding 
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Fig. 2. — All radial velocity measurements from Table [JJ folded 
on the orbital period: DA white dwarf uqa (open squares) and DQ 
white dwarf DjjQ (full triangles). The velocity amplitudes constrain 
the mass ratio to nearly unity (q = Mjjq/Meia = 1-1 ^ 0.1). 

SSO measurements to the set (N = 18) reduced the pe- 
riod uncertainty and allowed accurate phasing of the DQ 
velocity measurements as well (Table [2J with a period of 



P = 1.1697. 



-0.0057 
-0.0037 1 



Adding X-shooter data to the set (N = 19) resolved the 
solution into a triplet with each of the triplet members 
split into two possible periods. Focusing on the central 
peak, the first of the two solutions locates the X-shooter 
measurement near phase 0.15 (P — 1.1697d) while the 
second locates it symmetrically about the quadrature 
and near phase 0.35 (P = 1.1703d). Either solution 
is acceptable and the ambiguity does not affect the re- 
sulting orbital parameters for the DA and DQ stars. The 
same pattern is reproduced in the side peaks of the triplet 
but with a lesser significance. Table [2] lists the adopted 
orbital parameters. The initial epoch To corresponds to 
the lower conjunction of the DA white dwarf. Note that 
the period P and the initial epoch To are strongly anti- 
correlated. 

The period and velocity amplitudes constrain the to- 
tal mass of the system to (-Mda + Mdq) sin 3 i — 1.01 ± 
0.09 Mq, where i is the inclination of the orbital plane, or 
Mda + Mdq > 0.92 Mq. The combined velocity ampli- 
tude is Ada + Kdq = 202 km s" 1 . T he same amplitude 
estimated bv lVennes fe Kawkal ()2012f ). ATda+Adq > 320 
kms -1 , was in error with the correct constraint at only 
half of that value (> 160 kms -1 ). The consequences 
for the merger time-scale are discussed in the following 
Section. 

Figure [2] shows the radial velocity phased with the or- 
bital period. The X-shooter data are phased by adopting 
one of the two possible aliases (P = 1.1697 d). The dis- 
persion in Ha velocity measurements was only 10 kms -1 
and consistent with the estimated velocity errors. The 



Fig. 3. — Sets of co-added Ha spectra near quadratures. The Ho 
line core shifted by ~ 4A between the two sets. 



TABLE 2 
Binary parameters. 



Parameter 



Period 
Initial epoch (HJD) 
Semi-amplitude (DA) 
Mean velocity (DA) 
Semi-amplitude (DQ) 
Mass ratio 



Measurement 
7 +U.UUJiV 



-0.15 
+0.10 



P = 

T = 2455898.34 

K nA = 104.6 ± 3.7 'kms -1 
7DA = 69.3 ± 2.8 kms" 1 
Kdq = 97.4 ±5.1 kms" 1 

g = Af DQ /A/ D A = l.l±0.1 



DQ velocity curve is offset by —60 kms -1 relative to the 
DA velocity curve. The effect cannot be attributed to a 
difference in gravitational redshifts because the mass ra- 
tio is close to unity. Instead, the pressure shift of molec- 
ular band positions is a likely cause for this effect. In the 
following we utilize the DA mean-velocity in the calcula- 
tion of the systemic velocity and Galactic kinematics. 

Figure [3] shows two sets of co-added KPNO spectra, 
one near each quadrature. The velocity offset (Av « 180 
kms -1 ) proved readily measurable. 

If we adopt the spectroscopic mass estimate for the 
DA white dwarf based o n the measured (T e g,\ogg) from 
IVennes fe Kawkal (I2012I) . M DA = 0.829 ± 0.096 (fol- 
lowin g the mass-radius relations of lBenvenuto fc Althausl 
1999), and the orbital mass ratio (q = 1.1±0.1), then the 
total mass of the system is in the range 1.466-2.035 Mq. 
Consequently, the mass function constrains the inclina- 
tion to i = 52 — 70°, and the possibility of eclipses is 
ruled out. 

3.2. Age, evolution, and kinematics 

The cooling age of the DA white dwarf is estimated at 
7" C ooi = 1.9 ± 0.5 Gyr. A high mass for both white dwarfs 
i mplies relatively mas sive, Mj > 3 Mq (see a discussion 
in lFerrario et al.ll2005f). hence short-lived progenitors (~ 
0.48 Gvr. lSchaller et al.|[l992l ). constraining the total age 
of the system to less than 3 Gyr. 



4 



Vennes et al. 



Using our measured period (1.17 d) and mass ra- 
tio q = 1.1 ± 0.1, and constraining the DA mass to 
0.829±0. 096 Mp ), the m erger time is estimated at 10 11 yr 
following iRitterl (|198 6f). Therefore, the system is not a 
candidate type la supernova progenitor. 

The systemic velocity was obtained by subtracting 
the calculated gravitational redshift (7 gra v = 54 ± 6 
kms -1 ) from the mean velocity of the DA white dwarf 
(7da = 69 ± 3 kms" 1 ), 7 sys = 7da - 7g rav = 15 ± 7 
krns" 1 . Adopting a distance of 40 ± 6 pc, the systemic 
velocity and proper-motion measurements {fi a cos 8 = 
2 ± 6, us = -140 ± 6 mas yr" 1 , iSalim fc Gouldl 
120031) correspond to the Galactic velocity components 
(U, V, W) = (0 ± 8, - 2 7 ±6 , -6 ± 2) kms" 1 following 
Uohnson fe Soderbloml (119871). The syste m belongs to 
the thin-disk population ([Pauli et al.ll2006l ) in agreement 
with its young age (< 3 Gyr). 

4. SUMMARY AND DISCUSSION 

We have shown that the peculiar DQ white dwarf 
in NLTT 16249 is in a close 1.17 d orbit with a 
DA white dwarf compan ion. Close double - degen erate 
stars are common and Ma xted fc Marsh] (|1999() re- 
ported a fraction of 5 to 19% of close pairs in their 
radial velocity survey. White dwarfs with compos- 
ite spectra indicative of DA plus DB (Hei) or DA 
plus DC (cool He-rich) pairs are also relatively com- 
mon in large surveys (see, e.g.. iLimoges fc Bergeron! 
120101: iTremblav. Bergeron, fc Gia nninas 2 0111) . However, 
only two DA+DQ pairs are known: NLTT 16249 and 
SDSS J153210.04+ 135615.0 (NLTT 40489) wh ich was re- 
cently identified by iGiammichele et all (|2012[ ) . 

The velocity amplitudes in NLTT 16249 imply that 
the DQ is slightly more massive than its companion in 
which case it is more likely that the DQ formed first, both 
from ~ 3 Mq progenitors, unless stable mass transfer 
reversed the initial m ass ratio (for detailed scenarios see 
lNelemans~ et al. 200l|). The relatively long orbital period 
reported here rules out the prospect of a merger within 
a Hubble time. 

The DQ white dwarf is peculiar and the nitrogen 
concentration in it s atmosphere is, so far, unique. 
lAlthaus et~aH (|2005f ) followed the evolution of a 2.7M 
main-sequence star past the PG1159 spectroscopic stage 
and down the cooling track up to the DQ stage (0.6M©). 
The star follows a "born-again" He-flash loop that re- 
sets its evolution onto the post-asymptotic giant-branch 



(post-AGB), considerably modifying the chemical struc- 
ture of the outer layers that an ot herwise norma l post- 
AGB would possess. In particular, lAlthaus et al.l (2005) 
predicted measurable concentrations of nitrogen and oxy- 
gen in the atmosphere of DQ descendents of born-again 
stars. 

The nitrogen abundance measured in the DQ compo- 
nent of NLTT 16249 is qualitatively similar to the abun- 
dance predicted for the 10,500 K model, the lowest tem- 
perature considered bv lAlthaus et al.l (|2005| ). However, 
the C/N ratio for that particular model is close to 300, i.e, 
much larger than measured in NLTT 16249 (C/N rj 50). 
In the born-again context, the difference could be at- 
tributed to a different progenitor mass which could affect 
the outcome of the chemical evo lution, although such 
models are not available to us. lAlthaus et al.l ()2005l ) 
also predicted a concentration of oxygen between that 
of carbon and nitrogen (C/O ~ 100). It is possible to 
constrain the C/O abundance ratio in DQ white dwarfs 
using the near and far ultraviolet CO "Fourth-Positive" 
(4P) bands along with CN violet and C2 Swan bands. 
Finally, the 13 C/ 12 C isotopic ratio expected for the same 
evo lutionary track is well be low detection limits reported 
by iVennes fc Kawkal ([2012) in the case of the DQ in 
NLTT 16249. 

What would be the role, if any, played by binary in- 
teractions in shaping the chemical structure of the DQ 
progenitor in NLTT 16249 ? The DQ white dwarfs 
are not known to show nitrogen i n their spectra (see 
iDufour. Bergeron, fc Fontaine! 120 05V Therefore, it may 
not be a coincidence that the only known detection of ni- 
trogen in a DQ white dwarf is also that of a white dwarf 
in a close double-degenerate system. In that context as 
well, it remains to be shown that detectable traces of ni- 
trogen would be the expected outcome of mass-transfer 
or early envelope ejection. 
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